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ABSTRACT 

We present high-precision photometric light curves of five O-type stars observed with the 
refurbished Kepler satellite during its Campaign 0. For one of the stars, we also assembled 
high-resolution ground-based spectroscopy with the hermes spectrograph attached to the 1 .2- 
m Mercator telescope. The stars EPIC 202060097 (09.5V) and EPIC 202060098 (07V) ex¬ 
hibit monoperiodic variability due to rotational modulation with an amplitude of 5.6mmag 
and 9.3mmag and a rotation period of 2.63 d and 5.03 d, respectively. EPIC 202060091 
(09V) and EPIC 202060093 (09V:pe) reveal variability at low frequency but the cause is 
unclear. EPIC 202060092 (09V:p) is discovered to be a spectroscopic binary with at least 
one multiperiodic (3 Cep-type pulsator whose detected mode frequencies occur in the range 
[O.II,6.99]d"' and have amplitudes between 0.8 and 2.0mmag. Its pulsation spectrum is 
shown to be fully compatible with the ones predicted by core-hydrogen burning O-star mod¬ 
els. Despite the short duration of some 33 d and the limited data quality with a precision 
near lOO/tmag of these first K2 data, the diversity of possible causes for O-star variability 
already revealed from campaigns of similar duration by the MOST and CoRoT satellites is 
confirmed with Kepler. We provide an overview of O-star space photometry and give argu¬ 
ments why future K2 monitoring during Campaigns 11 and 13 at short cadence, accompanied 
by time-resolved high-precision high-resolution spectroscopy opens up the possibility of in- 
depth O-star seismology. 

Key words: Asteroseismology - Stars: massive - Stars: rotation - Stars: oscillations (includ¬ 
ing pulsations) - Techniques: photometry - Techniques: spectroscopy 
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1 INTRODUCTION 

Recent high-precision uninterrupted high-cadence space photom¬ 
etry implied a revolution in the observational evaluation of stel¬ 
lar structure models for various types of low-mass stars, covering 
almost their entire evolutionary path. Many results were obtained 
for sun-like stars and red giants undergoing solar-like oscillations 
excited stochastically in their convective envelope, thousands of 
which were monitored. The interpretation of their oscillation spec¬ 
trum is readily achieved by relying on existing methodology de¬ 
veloped for helioseismology, when extended with the interpreta¬ 
tion of gravity-dominated dipole mixed modes (e. g., iBedding et al.l 
I2OI itiMosser et alj2012l : ICliaDlin & Migliol2013h . Although many 
open questions r emain on the evolutio nary state of blue horizontal 
branch stars (e.g.. l0stensen et alj2012h . progress was also achieved 
from asteroseismology of a few subdwarf B stars and white dwarfs. 
These are evolved low-mass stars, the former of which have lost al¬ 
most their entire hydrogen envelope while passing through the he- 
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lium flash in a binary configuration (e.g., IVan Grootel et alJbOld : 
ICharpinet et alj 1201 ih . The latter are the compact remnants of 
low-mass stars that were very successfully studied from ground- 
based mm ag-precision asteros eismology more than two decades 
ago (e.g., ICorsico et al.l 1200 SL and references therein) but they 
were hardly monitored seismically from sp ace due to their faint¬ 
ness and fast oscillations of tens of seconds 1 0stensen et al.ll20ni : 
ISischotf-Kim & 0stensenl l201ll : iHermes et alj |201~ ). The heat- 
driven pressure and/or gravity modes detected in hundreds of 
intermediate-mass stars are much harder to analyse and interpret 
theoretically than the stochastically-excited solar -like pulsators, but 
major progress is also underway for suc h objects dKurtz et al.l2014l : 
ISaio et ^ .llZOlShlVan Reeth et alj2oi^ . 

From the point of view of chemical evolution of our Milky 
Way, supernova progenitors, i.e., stars with masses above ~ 9 Mq 
having an extended convective core and a ra diative envelope at 
birth, are the objects that matter l lMaedei|[2009l) . Despite their large 
importance for massive star evolution theory, O stars have hardly 
been monitored in high-cadence uninterrupted space photometry. 
If available, such monitoring was so far restricted to mmag preci¬ 
sion obtained by the WIRE, MOST, and Coriolis (SMEI instru¬ 
ment) satellites, while CoRoT observed six O stars during 34 d 
leading to a level of brightness variations with ~100^mag in pre¬ 
cision. Except for the SMEI light curves, which are seven to eight 
months long, all of the available time series have a duration of the 
order of a month only. 

Unfortunately, O stars were absent in the nominal Kepler 
Eield-of-View (FoV), hence we do not have long-term /rmag-level 
precision photometry of O stars at hand, preventing detailed aster- 
oseismic calibrations of massi ve star mode ls as it was achieved for 
the less massive B stars te.g. lAert3l2015L for a summary). Sim¬ 
ilarly to the O stars, B stars have an extensive convective core 
but they do not exhibit a strong radiation-driven wind. The two 
B stars that were seismically modelled based upon period spac- 
ings of their high-order dipole gravity-mode oscillations are ultra- 
slow rotators (vsini < lOkms^'). Their modelling required the 
inclusion of extra ditfusive mixing, in addition to core overshoot¬ 
ing, to bring the theoretical models in agreement with the yumag- 
precision s eismic data of durati on five months (for the B3V star 
FID 50230, iDegroote et all 20 l^ and four years (for the B8V star 
KIC 10526294, Moravveii et alJl2oT5h . It would be highly benefi¬ 
cial to perform similar seismic inference studies for O-type stars, 
whose theoretical models are the most uncertain of all the mass 
ranges, while they are mainly responsible for the global chemical 
enrichment of the Milky Way. 

Two rapidly rotating O stars were monitored by the 
MOST satellite. It concerns fPer (OTIIInf) with vsini 
200kms“* revealing variability due to rotational modulation 
with a rotation period of 4.18 d, in the absence of oscillations 
jRamiaramanantsoa et aljl2014h and /(Oph (09.5Vnn, vsini 
400 km s ') whose WIRE, MOST, SMEI photometry led to the 
detection of multiperiodic variability due to oscillations with fre¬ 
quencies in the range [1,10] d“' and amplitudes up to =: 10 mmag. 
These modes turn ou t to be variable on a time scale of some 100 d 
dHowarth et alJl2Ql4ll . The rapid rotator fPup was monitored by 
SMEI during several long-term runs, but these data revealed only 
two dominant frequencies, whose cause and identification remain 
unclear; they might be con nected with non-adiabatic gravity modes 
dHowarth & Stevenil2014h . The variability derived from the space 
photometry of these three rapid rotators is in line with assess¬ 
ments from their high-precision time-resolved line-profile variabil¬ 
ity. It is also in agreement with high-precision time-resolved spec¬ 


troscopy of the rapidly rotating 09Vp st ar HD 93521, whose multi- 
periodicity has a yet unc lear origin (e.g.. lFullerton. Gies. & BoltonI 
1 19961; iRauw 61^120081 vsini =! 400kms *). Seismic modelling 
of rapidly rotating O stars was not achieved so far. 

One of the CoRoT short runs (SRa02) was devoted to the 
monitoring of six O-type stars during 34 d in th e asteroseismol- 
ogy C CD of the mission, having a cadence of 32 s dAuvergne et al.l 
I2OO9I) . Among these six stars were two binaries. Plaskett’s star 
(HD 47129) is a high- mass interacting bin ary (07.5I-I-06II) with 
a magnetic secondary dOrunhut et al.ll20f^ . Its CoRoT light curve 
revealed hints of gravity-mode oscillations with frequencies near 
~0.8d“‘ and harmonics, as well as rotation al modulation with 
multipl es of the orbital pe riod of 14.39625 d dUinder et al.ll2008l) 
— see iMahv et^ d201 ih for a discussion. A seismic interpre¬ 
tation was not possible for this complex system. The eccentric 
(e = 0.59 ± 0.02) long-orbit (Porbit = 829 ± 4d) binary HD 46149 
revealed a primary with stochastic pressure-mode oscillations with 
frequencies in the range 3.0 to 7.2 d“'. These correspond to a reg¬ 
ular frequency pattern of separation 0.48 d“‘. Unfortunately, these 
p modes did not l ead to an unambiguo us interpretation in terms 
of seismic models |Degroote_et^ 20 id). Fur ther, the “red noise” 
power excess found by Blomme et alj d201 ih in the CoRoT am¬ 
plitude spectra of the three moderately rotating O stars HD 46150 
(05.5f, vsini 100kms-‘), HD46223 (05f, vsini == 100kms-'), 
and HD46966 (08.5V, vsini 50kms“') were recently in- 
terpreted in terms of convectively-driven internal gravity waves 
dAerts & Rogersll2015l) . 

The only seismic inference was achieved for the 09V star 
HD 46202, which revealed /t Cep-type pressure-mod e oscillations 
with frequencies in the range from 0.5 to 4.9 d“‘ dBriguet et al.l 
l201lh . Seismic modelling of this star led to a mass M = 24.1 ± 
0.8 Mo, an age of 4.3 ± 0.5 Myr, and a core overshooting parameter 
of 0.10 ± 0.05 expressed in units of the local pressure scale height 
and using a step-function formulation. 

To remedy the la ck of seismic calibr ation of massive star mod¬ 
els, we defined a K2 dHowell et alj|20l3) observing programme of 
O-type stars with the aim to monitor a sample of such objects at 
ytimag precision during 3-month runs. The present work reports on 
the first results of this programme and otfers a proof-of-concept of 
K2’s capabilities for future O-star asteroseismology. 


2 EXTRACTION OF THE LIGHT CURVES 

We proposed the five O-type stars listed in Table[T]for photometric 
monitoring with K2 in Campaign 0. We developed a photometric 
reduction method that accounts for the spacecraft (re)pointing and 
the non-uniform pixel response of the K2 CCDs, starting from the 
pixel frames in the MAST archive (release of 7 Nov. 2014). The 
time series of Campaign 0 contain two large gaps related to safety 
events of the spacecraft. We only used data after the largest safety 
event, i.e., after the Kepler Barycentric Julian Day 1939.1, provid¬ 
ing us with photometric measurements covering 33.038 d. These 
data have a frequency resolution corresponding to the Rayleigh 
limit of 1/33.038 d= 0.0302 d“* (0.35;uHz), which represents an 
upper li mit for the freque ncy uncertainty of any periodic signal in 
the data dAerts et al.l2010l . Chapter 5). Further, we ignored all mea¬ 
surements taken during a thruster firing or a repointing event of the 
spacecraft. 

With the aim to exclude as much as possible instrumental ef¬ 
fects in the light curves, we first defined an optimal constant mask 
customized to each of the five targets from visual inspection, such 
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Figure 1. K2 light curves (top) and their amplitude spectra (bottom) of the five m onitored O stars. T he black circles and thin black lines are derived from our 
reduction of the raw pixel data, while the grey crosses and thick grey lines are the lVanderburd i2014l) light curves. 
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Table 1. The five O stars observed in Campaign 0 of K2, with some information from the online SIMBAD database and the EPIC catalogue. The dominant 
frequency in the K2 data is listed, along with its amplitude. 


EPIC ID 

Other ID 

RA (J2000) 

DEC (J2000) 

V (mag) 

SpT 

Frequency (d 

Amplitude (ppt) 

202060091 

HD 44597 

06 23 28.538 

-t20 23 31.63 

9.05 

09V 

— 

— 

202060092 

HD 256035 

06 22 58.243 

+22 51 46.15 

9.21 

09V:p 

3.4248 

1.827 

202060093 

HD 255055 

06 19 41.647 

+23 17 20.22 

9.40 

09V:p(e) 

— 

— 

202060097 

Cl* NGC 2175 H 98 

06 08 31.999 

+20 39 24.01 

13.7 

09.5V 

0.3801 

5.209 

202060098 

2MASS J06245986-t2649194 

06 24 59.870 

+26 49 19.42 

15.11 

07V 

0.1988 

8.594 


that blending with nearby targets is avoided as much as possible. 
We show these masks in the left panels of Figs lAll to IA5I The 
counts for each pixel in these masks were summed up to provide 
the time series of the raw photometric measurements. To account 
for the varying background during the measurements, we used it¬ 
erative sigma clipping. All pixels of the frame at a given time step 
were deselected when their signal was larger than 1.2 times the me¬ 
dian signal of the frame. This procedure was iterated until we found 
a stable median for that time step and no further deselection was 
needed. The achieved median was then taken to be the background 
signal in a given pixel in the frame. We subtracted this background 
level from the raw light curve, accounting for the considered num¬ 
ber of pixels in the mask. 

Outliers were identified and removed by self-flat-fielding 
through spline fitting, along with further detrending of the pho¬ 
tometry by third-order polynomial fitting, in an iterative scheme. 
In each step, the photometry was divided by the polynomial fit, 
which allowed us to work with relative brightness units instead of 
electrons per second. Five iteration steps turned out to be sufficient. 

This raw light curve was then used as a starting point to ac¬ 
count for the instrumental effects caused by K2’s attitude control 
system, which involves a roll manoeuvre every six hours. Due to 
that, we had to redetermine the position of the star on the CCD 
and connect it to the instrumental brightness variations that oc¬ 
curred at that position. In order to derive an appropriate relation 
for the correction factor versus star position, we determined the po¬ 
sition of the star on the CCD by fitting a 2D Gaussian function 
to t he photometry, representing th e point-spread-function. Follow¬ 
ing |^ndeAurg|&7ohns^ ( l2014h . we rotated that 2D position ac¬ 
cording to its largest eigenvalue and fitted a fifth-order polynomial, 
allowing us to define an arclength connected with the instrumental 
brightness for each time step. We then assumed that all the power in 
the Fourier transform of the time series of the measured brightness 
centred at 4.08 d“* and its higher harmonics at 8.16, 12.24, 16.32, 
and 20.40 d“* had instrumental origin. A bandpass filter in the 
Fourier domain was constructed using the bands 7X4.08 d“' ±5/^5, 
with /res the Rayleigh limit and j = 1,...,5. All the power in 
these bands was assumed to be due to instrumental effects, while 
all power outside these bands was considered as stellar signal. The 
inverse Fourier transform then allowed to connect the instrumental 
brightness variations to the position of the star on the CCD. This 
method worked well for four of the five stars, EPIC 202060093 suf¬ 
fering from a too small onboard mask and saturation of many pixels 
leading to leakage outside of the frame. The light curves are shown 
with black circles in the upper panel of Fig.[T] These flux variations 
are expressed in parts-per-thousand (ppt)[j 

In order to compare our reduction scheme for the individual 


stars and its effect iveness, we also us ed the reduced light curves as 
made available bv IVanderburd ( l2014h . release date: 10 April 20150 
This is an updated reduction based on their self-flatfielding ap¬ 
proach combined with high-pass filtering, designed from the view- 
point of optimal exoplanet h unting, as previously developed by 
IVanderburg & Johnsonl ( l2014h to which we refer for details. While 
our mask determination and noise treatment was done manually 
and tuned to each star individually and locally, IVanderburd ( l2014h 
used a semi-automated mask determination and made use of the in¬ 
strumental effects of thousands of stars. Our method is tuned to get 
the optimal signal, while their method is focused on reducing the 
global noise as much as possible. The masks used by Vanderburg 
are shown in the right panels of Figs lAll to lASI 

Despite the quite different masking, the two independent re¬ 
duction methods give consistent results in terms of frequencies, the 
major difference occurring in the lowest frequency regime of the 
Fourier transform. This is similar to the difference already encoun¬ 
tered for light curves of the nominal Kepler mission when reduced 
globally from the viewpoint of optimal planet hunting versus in¬ 
dividually and tuned towards optimal asteroseismic applications. 
For the latter, it is typically advantageous to consider more pixels 
than done in the standard exoplanet pipeline, because the aim is to 
gather as much signal as possible for a whole range of frequencies, 
given the multiperiodic character of the variability. For exoplanet 
hunting, however, one aims to find singly periodic low-amplitude 
signal and therefore avoids too noisy pixels. Masking for astero- 
seismology as we have done here avoids complicated detrending 
and filters out long-term i nstrumental frequencies more easily, as 
illustrated and discussed in iTkachenko et all l l2013ll . but is not op¬ 
timal in terms of global noise properties. 

A comparison between the light curves resulting from the two 
independent reduction methods for the five O stars considered here 
is given in Fig.[T] In general, the agreement is good, except for 
EPIC 202060098 and EPIC 202060093. For the former star, the 


dominant frequency is fully consistent but its amplitude differs in 
value, due to the different shape and position of the chosen masks 
in the two methods. With our mask, we have avoided the additional 
bright target contributing to the overall light curve. This explains 
the higher amplitude compared to the “Vanderburg” light curve, 
while the latter leads to less power at low frequencies (cf. bottom 
panel of Fig.[T](. For EPIC 202060093, both light curves suffer from 
the saturation and leakage across the CCD, but again the frequency 
amplitude is higher when derived from our light curve. In the rest of 
the paper, we used our reduced light curves to derive the frequen¬ 
cies and amplitudes of the five target stars, after having verified that 
all the listed frequencies are consistently recovered from those two 
versions. 


* The conversion factor between brightness variations expressed in mmag 
and flux variations expressed in ppt amounts to 2.5 logjo e = 1.0857 


https://WWW.cfa.harvard.edu/~avanderb/ 
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Table 2. The frequencies and amplitudes of EPIC 202060092 determined 
from Fourier analysis followed by non-linear least-squares fitting. The com¬ 
putation of the errors is explained in the text. The signal-to-noise ratio 
(SNR) was computed as the amplitude divided by the noise level of the 
residual light curve averaged over [0,10] d“*. 


ID 

Frequency (d ^) 

Amplitude (ppt) 

SNR 

Remark 

/l 

3.4243 

± 

0.0054 

1.827 

± 

0.079 

16.9 


/2 

3.6663 

± 

0.0066 

1.352 

± 

0.084 

12.5 


h 

6.9853 

± 

0.0088 

1.091 

± 

0.104 

10.1 


/4 

3.3173 

± 

0.0113 

0.938 

± 

0.103 

8.7 


fs 

3.5583 

± 

0.0089 

0.998 

± 

0.105 

9.2 


fi 

1.7453 

± 

0.0092 

0.941 

± 

0.083 

8.7 

/3/4 

fl 

1.6693 

± 

0.0085 

0.953 

± 

0.106 

8.8 

/4/2 

h 

1.7193 

± 

0.0091 

0.784 

± 

0.106 

7.3 

/i/2 

/9 

3.4863 

± 

0.0113 

0.646 

± 

0.096 

6.0 

fi/2 

/lO 

0.2223 

± 

0.0113 

0.705 

± 

0.100 

6.5 

2fn 

/ll 

0.3313 

± 

0.0125 

0.585 

± 

0.104 

5.4 

3fn 

/l2 

0.1123 

± 

0.0129 

0.566 

± 

0.074 

5.2 

fl- h = h- fs 

/l3 

1.3223 

± 

0.0134 

0.535 

± 

0.095 

5.0 



3 FREQUENCY ANALYSIS OF K2 DATA 

We computed the Lomb-Scargle periodograms of the two versions 
of the light curves. These are shown in the lower panel of Fig.[T] 
The thin black lines are for our reduced light curves while the thick 
grey lines are for the lVanderburd (1201 4) light curves. 

For the stars EPIC 202060097 and 202060098, we find one 
isolated significant frequency corresponding to a period of 2.631 d 
and 5.030 d, respectively (Table[T). We interpret their variability as 
due to rotational modulation. 

The cause of the variability of EPIC 202060091 and 
202060093 is less clear. Despite the limitations of the masks for 
EPIC 202060093 (Fig IA3b . its frequency spectrum is similar to the 
one of EPIC 202060091, but there is a factor more than ten dif¬ 
ference in the level of variability. None of these two stars shows 
clear periodic variability with isolated frequencies. Nevertheless, 
they display several low-frequency peaks that stand out of the noise 
level without being formally significant. This is somewhat similar 
to the frequency spectra of three O stars observed with the CoRoT 
satellite, revealing red noise power excess at low frequency that was 
recentl y interpreted in terms o f convectively driven internal gravity 
waves ( lAerts & Rogersl2015h . However, we currently consider this 
interpretation as speculative. 

EPIC 202060092 turns out to be a multiperiodic pulsator com- 

S iatible with heat-d riven oscillation modes, similar to HD 46202 
Briquet et al ■1201 ih . We performed iterative prewhitening and find 
13 significant frequencies when we adopt the conservative crite¬ 
rion of considering a frequency to be significant when its amplitude 
reaches above five times the noise level in the residual light curve 
averaged over the frequency range [0,10] d^F The frequencies and 
amplitudes along with their errors are listed in Table[2]while the fit 
to the light curve based on those 13 frequencies is shown in Pig.|2] 
Eor the computation of the frequency and amplitude errors, we took 
into account that the formal errors resulting from a non-linear least 
squares fit are only in agreement with those obtained in the Eourier 
domain in the case of uncorrelated data w ith white noise an d with 
sufficiently high frequency resolution (e.g., lAerts et alj20ld . Chap¬ 
ter 5). Here, we encounter two complications: the data have only a 
short time base of some 33 d, le ading to a limited resolving power 
of 1.5/33.038 d = 0.04541 d^' dUoumos & Deeminglfl^ . Eur- 
ther, as is usually the case for highly sampled space photometry, the 
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Figure 3. Selected H and He lines in the highest SNR HERMES spectrum 
of EPIC 202060092. The blue line is the prediction for a NLTE atmosphere 
model with parameters listed in the text. 


data may be correlated. This requires incorporation of a correction 
factor in the error estimates of the frequenc ies and their amplitude s 
dSchwarzenberg-Czemvl l2003h . Eollowing IPegroote et all d2009h . 
we computed this correction factor to be 2.51 and obtained the er¬ 
ror estimates listed in Table|^ The fit shown in Eig.l^is generally 
good, but not perfect. This is due to there being additional frequen¬ 
cies whose amplitudes are between 3 and 5 times the noise level. 
Due to the limited frequency resolution, we are unable to pinpoint 
their value with respect to the frequencies already listed in Table|2] 
and this explains why a level of unresolved beating still occurs in 
some parts of the light curve (bottom panel of Eig.[2]l. 


4 FOLLOW-UP STUDY OL EPIC 202060092 

The nature of the frequencies /e to /b listed in Tablel^is unclear. 
Given the limited frequency precision, they could either be identi¬ 
fied as sub-multiples or combinations based on /, to f^, or be due 
to independent g modes occurring in the densely-packed gravity¬ 
mode frequency regime, or both. The interpretation of all the de¬ 
tected frequencies in terms of identification of their degree radial 
order «, and azimuthal order m requires additional observational 
information, such as the rotational frequency of the star. We also 
note that the light curve of EPIC 202060092 shows a marked dip 
near day 1966 and remarkable beating from days 1946 to 1948 and 
to a lesser extent also from days 1963 to 1965 (Eig.|2l where the 
listed and shown dates are with respect to Kepler Barycentric Ju¬ 
lian Data, indicated here as BID2454833.0). With the aim to test 
if the dips could be connected with binarity and to derive the pro¬ 
jected rotation velocity of the pulsator, we gathered spectroscopic 
measurements. 


4.1 EPIC 202060092 is a spectroscopic binary 

The 1.2-m Mercator telescope is dedicated to the long-term mon- 
itoring of variable stars, including hea t-driven pulsators (e.g., 
IPe Cat et al.l 200j|_ Cuy perset alj l2009h and evolved low-mass 
stars (e.g., Van Winckel et al.l 20091) . A large fraction of the tele¬ 
scope time in the recent years has been used for the spectroscopic 
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Figure 2. Fit (grey line) to the K2 light curve of EPIC 202060092 (black circles repeated from the central panel of Fig.0 based on the 13 frequencies listed in 
Table[2] 
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Figure 4. Phased radial velocity measurements of EPIC 202060092 de¬ 
duced from its Balmer lines. The residual K2 light curve after prewhitening 
the fit in Fig.[2is indicated in the top part in ppt according to the same phase 
and shifted up with value 15 for visibility reasons. 


monitoring of bright Kepler targets (e. 

g.. Lehmann et al. 

2011; 

Tkachenko et al.l 20121; iReck et alJl2014 

; iNiemezura et al. 

2Q15; 

Van Reeth et akluOlsb. In view of its photometric behaviour and of 


recent theoretical predictions that m ore than 70% of all 0-type stars 
occur in binaries Isana et d ]l2012h . we added EPIC 202060092 to 
the obser ving programme fo r measurements with the hermes spec¬ 
trograph (iRaskin et al.ll201lh . 

A total of 22 high-resolution (HRF mode, resolution of 
85 000) spectra were taken so far, with integration times between 
600 s and 1800 s between November 2014 and March 2015. The 
signal-to-noise ratio (SNR) per exposure ranged from 30 to 60. The 


raw exposures were treated with the HERMES pipeline, including 
cosmic-hit removal, merging of orders, and barycentric correction. 
Subseq uently, all spectra w ere rectified following the method out¬ 
lined in lPapics et alj ( l2013h . 

Comparison with rotationally broadened synthetic spectra 
derived from plane-pa rallel NLTE 0-star atmosphere models 
jLanz & Huben^ l200l OSTAR2002 grid) led to the estimates 
Tcff - 35 000 K, log g 4.5 dex, vsin i ^ 270kms“* for solar 
metallicity and a fixed microturbulent velocity of lOkms^*. The 
agreement between the theoretical line spectrum and the measured 
profiles for the highest SNR spectrum is illustrated in Fig.|3]for a 
few selected lines. Despite the modest SNR, it can be seen that 
structure occurs in some of the observed He lines. This may be 
connected to the pulsations but higher SNR is needed to firmly es¬ 
tablish this. 

Given the scarcity of spectral lines due to the large broaden¬ 
ing and limited SNR, we estimated the radial velocity values from 
fits to the detected H and He lines. The outcome is very similar for 
the Her, H/3, Hy, and He lines, so we averaged their radial velocity 
estimates and computed the standard deviation as a proxy for the er¬ 
rors. We detect radial velocity variations with a peak-to-peak value 
around 30 km s“'. Given the low amplitude of the photometric vari¬ 
ations (cf. Table|2]l we do not assign these radial-velocity variations 
to pulsations, but rather interpret them in terms of binarity. Glob¬ 
ally, the velocities from the Balmer lines are consistent with those 
based on the dominant He lines, such as HeI4922A, HeI5875A, 
Hel 6678A, Hell 4541 A, and Hell 5411 A, but for some of those we 
encounter too large uncertainties to achieve conclusive values, as 
already illustrated by Fig.[3] This might point towards a contribu¬ 
tion of more than one star to some of these lines, but this needs 
further study from new spectra with higher SNR. 

Our current spectroscopy does not allow to deduce a high- 
precision orbital solution yet. Nevertheless, we clearly establish 
spectroscopic binarity and find a rough preliminary period estimate 
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of ^orb ~ 50d from the Balmer lines. These data were folded with 
the most likely period in Fig.|4l where we also show the phased 
residual K2 light curve after prewhitening the 13 frequencies listed 
in Table[2 (i.e., the black circles minus the grey line in Fig.^. 
We conclude that, with the current estimate of the orbital period 
of the binary, we cannot associate all the dips in the K2 photom¬ 
etry with eclipses and rather ascribe them to unresolved beating 
of pulsation modes, except perhaps the dip near day 1966. Large 
beating patterns as the ones seen near days 1946, 1955, and 1964 
(Fig.|3 have indeed been detected before in multipe riodic massive 
pulsa t ors, b oth in high-cadence ground- based (e.g. Handle r et al.l 
|2004l200^ and space photometry fe.g.. |PaDics et alj|2012tl . With 
the long-cadence sampling at 29 min, the unresolved nature of the 
light curve is not so surprising, because the fastest pulsation mode 
is sampled only seven times per cycle. While such sampling works 
fine for the analysis and interpretation of long-duration light curves 
of months to years, there is a limitation to resolve beating patterns 
for a data set of only 33 d. The beating between frequencies /i 
and / 4 , as well as the one between fn and /s, occurs with a period 
of 9.1 d and the beating phenomena we measure are separated by 
roughly this value, occurring near days 1946, 1955, and 1964 (top 
panel of Fig.|3- The nature of the isolated dip in the light curve 
at day 1966 remains unclear. That decrease in brightness is not 
accompanied by local increases as one would expect from multi¬ 
periodic beating among linear pulsation modes and could therefore 
still be caused by binarity. Indeed, according to Fig.|4] the dip oc¬ 
curs near orbital phase ~ 0.5 and, given the large uncertainty on 
the current estimate of the orbital period, could be caused by an 
eclipse. Unfortunately, there are no FIipparcos data data available 
to add to the K2 photometry for this star. 


4.2 First modelling attempts 

The five dominant frequencies listed in Table|2]have values typ¬ 
ical of pressure modes in an 09 mai n-sequence star as alr eady 
found for FID 46149 and FID 46202 bv IPegroote et al.l ( 1201 Oh and 
iBriguet et alj ( 1201 ih . respectively. Part of the theoretical radial and 
zonal dipole or quadrupole p mode spectrum for 0-type stars is 
shown in Fig.j^for a mass of 35 Mq (left panel) and 30 Mq (right 
panel). This figure was constr ucted from the grid of stellar and pul- 
sational models descri bed in DegrootejetalJ (l2010h . whose input 
physics is discussed in lBriguet et alj ( 1201 ll) . to which we refer for 
details. The two stellar models whose zonal p modes are shown in 
Fig.|5]have an initial hydrogen fraction of Xc = 0.715 and a step- 
function core overshoot parameter of o-ov = 0.2, expressed in units 
of the local pressure scale height. As we have shown in the previ¬ 
ous Section, EPIC 202060092 is situated between the two vertical 
dashed lines. These correspond to ages of 3.53 and 3.51 Myr (hot¬ 
ter temperature limit) and 4.39 and 4.86 Myr (coolest temperature 
limit), for 35 Mq and 30 Mq, respectively. 

Even though relatively large uncertainties occur for the atmo¬ 
spheric parameters of the star, the measurement of v sin i is fairly 
robust. It immediately excludes fn = 0.1123d“' to be the rota¬ 
tional frequency. Indeed, the models with 35 Mq and 30 Mg whose 
low-degree p modes were shown in Eig.j^have radii of 15 Rg and 
12.5 Rg, respectively, for the measured Tat- We find lower limits for 
the rotational splitting within p mode multiplets of roughly 0.36 d“* 
and 0.43 d“' for 35 Mg and 30 Mg, respectively. 

Even though the white-light CoRoT photometry of HD 46202 
by itself did not lead to mode identification in the absence of multi¬ 
plets in the frequenc y spectrum, the same situation as we have here 
for EPIC 202060092. iBriouet et alj ( l201 ih managed to perform for¬ 


ward seismic modelling from frequency matching and pinpointed 
the four free parameters (M, X, Z, ttov) of HD 46202 with high pre¬ 
cision, for the input physics that went into their dense grid of the¬ 
oretical models. Their frequency precision was much higher, typ¬ 
ically between lO^”* and 10“^ d“‘, than what we can deduce from 
the current long-cadence K2 data. This major accomplishment for 
HD 46202 was owed to the CoRoT sampling rate of 32 s, which 
delivered a light curve of more than 80 000 data points rather than 
only 1 600 for EPIC 202060092 from K2. Despite similar duration 
of the monitoring, this implies a factor ~ VSO better frequency pre¬ 
cision. Moreover, their study was based on twice as many spectra of 
a single star of one magnitude brighter than EPIC 202060092, im¬ 
plying a co-added averaged spectrum of SNR above 500. Einally, 
the ten times lower v sin i = 25 ±1 km s“' of HD 46202 compared 
to EPIC 202060092 also implied modest rotational splitting values 
such that forward modelling assuming zonal modes made sense. 

Eor the case of EPIC 202060092, which has a similar oscil¬ 
lation spectrum and hence the same potential, we encounter the 
limitations of the poor frequency precision and of not knowing the 
azimuthal order m of the dominant modes because the fast rotation 
introduces frequency shifts of at least 0.3 d“' when m 0. It can 
be seen from Fig.l^that, even for only two masses and one set of 
values (A,Z, ttov), several of these models fulfill the requirements 
of the five measured frequencies, keeping in mind that large yet un¬ 
known frequency shifts due to rotational splitting must be allowed 
for. A similar situation of unknown rotat ional shifts was enc oun- 
tered in the NGC 884 cluster modelling bv ISaesen et al.l (l2013h . but 
in that study some seismic constraints could still be achieved thanks 
to the demand of an equal age and metallicity of the stars in the 
cluster. 

With the present K2 photometry and HERMES spectroscopy, 
forward seismic modelling is not yet possible for EPIC 202060092. 
We do note that even the identification of m) of just one of its 
modes would break most of the degeneracy in the seismic grid. As 
an example, identification of fn as a radial mode would imply it is 
the fundamental, while is then the third (for 35 Mg, left panel of 
Fig.|5) or fourth (for 30 Mg, right panel of Fig.[3 radial overtone. 
At present, however, such inference remains speculative. Neverthe¬ 
less, the models imply that, if /| to /j are due to low-degree low- 
order modes, then the mass of the star cannot be far below 30 Mg. 

Eirm conclusions on the properties of the orbit of 
EPIC 202060092 and on the nature of its companion require addi¬ 
tional and long-term time-resolved high SNR spectroscopy. Given 
its potential for future seismic modelling, we plan to acquire such 
data with the aim to pinpoint the orbit and characterise the two 
components with high precision, as well as to detect and identify 
the dominant oscillation mode(s). This will require spectral disen¬ 
tangling and the a pplication of sophistica ted mode identification 
methodology (e.g.. iTkachenko et al.ll2014L for an example of an¬ 
other potentially interesting pulsating massive binary). 


5 STATE OF AFFAIRS OF HIGH-PRECISION O-STAR 
SPACE PHOTOMETRY 

The five O-type stars monitored by K2 so far confirm earlier find¬ 
ings based on high-precision uninterrupted space photometry, i.e., 
the diversity of the variable character of these stars is large. Table[3] 
provides an overview of all fourteen O-type stars that have been 
monitored with SMEI, MOST, CoRoT, and K2 and their dominant 
cause of the photometric variability. For several of the stars, vari¬ 
ous causes of variability act simultaneously. The secondary of the 
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Figure 5. Adiabatic pulsation frequencies of radial (blue circles), dipole (black crosses), and quadrupole (red triangles) zonal modes for a stellar model with 
M = 35 M© (left) and with M = SOM© (right). EPIC202060092 is situated in between the two vertical dashed lines, which indicate the spectroscopic Icr-range 
of the measured Jeff. The top jc-axes indicate the central hydrogen fraction, where the value at birth was Xc = 0.715. The horizontal grey lines indicate the 
five dominant independent frequencies, where the line thickness represents the frequency precision. 


Table 3. Overview of 0-star variability deduc ed from high-cadenc e high-precision space photometry with the dominant causes of the variability listed, where 
IGW stands for “Internal Gravity Waves” fe.g. lRogers et alj2013h . The duration of the longest uninterrupted light curve is listed. 


HD or EPIC 

V 

mag 

SpT 

Data Source 

Duration 
ofLC (d) 

vsin i 
kms“^ 

Variability 

Cause(s) 

Reference to 

Space Photometry 

HD 46223 

7.28 

04V(f) 

CoRoT 

34.3 

100 

IGW 

Blomme et al. (201 11 

HD 66811 

2.24 

04I(n)fp 

SMEI 

236 

400 

non-adiabatic g“ mode? 

Howarth & Stevens (20141 

HD 46150 

6.73 

05V(f) 

CoRoT 

34.3 

80 

IGW 

Blomme et al. (201 11 

HD 24912 

4.06 

07IIInf 

MOST 

30.0 

200 

Rotational Modulation 

Ramiaramanantsoa et al. ('2014) 

EPIC 202060098 

15.1 

07V 

K2 

33.1 

? 

Rotational Modulation 

this work 

HD 46966 

6.87 

08V 

CoRoT 

34.3 

50 

IGW 

Blomme et al. (201 1) 

HD 47129 

6.06 

08 III/I+07.5III 

CoRoT 

34.3 

75 & 300 

Binarity, Rotation 
Magnetic secondary 

Mahvetal. 120111 

HD 46149 

7.61 

08.5V((f))+B? 

CoRoT 

34.3 

30 &? 

Rotational Modulation 
& Stochastic p modes 

Degroote et al. 120101 

HD 46202 

8.19 

09V(f) 

CoRoT 

34.3 

25 

Heat-driven p modes 

Briauet et al. (20111 

HD 44597 

9.05 

09V 

K2 

33.1 

~30 

IGW? 

this work 

HD 256035 

9.21 

09V:p+? 

K2 

33.1 

270 

Heat-driven modes 

this work 

HD 255055 

9.40 

09V:p(e) 

K2 

33.1 

~ 40 

IGW? 

this work 

HD 149757 

2.56 

09.5Vnn 

WIRE,MOST,SMEI 

232 

400 

Heat-driven p modes 

Howarth et al. (20141 

EPIC 202060097 

13.7 

09.5V 

K2 

33.1 

? 

Rotational Modulation 

this work 


binary HD 47129 is the only object in the list known to host a mag¬ 
netic field and it is variable. Nevertheless, rotational modulation 
connected with chemical and/or temperature spots occurs for 4 of 
the 14 stars and is usually due to some level of magnetic activity 
even though the fields may he too weak or too complex to give 
detectable longitudinal field components. 

From the viewpoint of seismic modelling with the aim to im¬ 
prove the input physics of massive star models, one needs to se¬ 
cure high-precision frequency values of several pulsation modes. 
While this is not impossible an d has been achieved for a rapidly- 
rotating pulsating Be star fe.g. iNeiner et alj|2012ll . it turned out 
hard to achieve so far in the case of i nternal gravity waves d riven 
by the convective core in O stars fe.g. lAerts & Rogeral2015h . The 


requirement to detect oscillation frequencies with high precision 
of ~ 0.0001 d“' can hest be achieved in the case of self-driven 
modes caused by an opacity mechanism active in the partial ion¬ 
isation layers of the iron-group elements, because they have long 
lifetime. HD 4620 2 is the only O-typ e star that has been modelled 
seismically so far teriguet et al.l201 ll) . While this led to the deriva¬ 
tion of its fundamental parameters with higher precision than any 
other method delivered so far, it was not possible to probe its in¬ 
terior structure properties at a level necessary to improve the input 
physics of the models, due to the too short time base of the space 
data. 

Heat-driven O-type pulsators hold the same potential of 
deep seismic sounding of their interior structure than was re- 
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cently achieved by iMoravveii et alj l l2015l) for the B8.3V star 
KIC 10526294. This 3.2 Mq star provided seismic evidence for dif¬ 
fusive mixing in the stellar envelope at a level of log Dmix between 
1.75 and 2.00 dex, in addition to core overshooting. A major im¬ 
provement for evolutionary models of massive stars requires the 
observational calibration of the overall mixing properties in their 
interior, because theoretical considerations lead to log Dnjx values 
differing by many orders o f magnitude and are therefore of limited 
value for the moment te.g.. lMathis. Palacios. & Zahnl2004}) . Aster- 
oseismology of O star pulsators is the best method to make progress 
for the tuning of models of the most massive stars in the Universe, 
but it requires long-term high-precision time-resolved photometric 
and spectroscopic monitoring. As we have shown in this work, the 
K2 mission has a major role to play here and EPIC 202060092 was 
an optimal target among five measured ones, but its complex beat¬ 
ing pattern was not sufficiently sampled in terms of cadence, preci¬ 
sion, and duration of the photometry achieved during Campaign 0. 
Moreover, it turned out to be a spectroscopic binary and we need 
to understand the contributions of each of its components to the de¬ 
tected variability. Requirements for its future seismic modelling are 
the accurate determination of its orbital properties and of the fun¬ 
damental parameters of its binary components, along with a higher 
precision of the pulsation frequencies. Given its potential, spectro¬ 
scopic monitoring will be continued in the coming years with the 
aim to unravel the orbital motion and to attempt detection of its 

oscillation mo des in long-term residu al s pectroscopy. _ 

Following iDegroote et al.l ( l2010ll and lBriguet et al.l ( l201lh . our 
current work is a successful proof-of-concept study to perform fu¬ 
ture 0-star seismology. This is possible with the optimised K2 mis¬ 
sion in combination with ground-based spectroscopy. Indeed, as of 
K2’s Campaign 3, the improvement on its performance in terms of 
pointing and precision of the photometry, as well as in the dura¬ 
tion of the campaigns (~ 80 d, i.e., twice to three times as long 
as the current study) holds the potential to deliver the appropriate 
photometric data needed for in-depth modelling of 0-star pulsa¬ 
tions, whatever their excitation mechanism. A light curve at short 
cadence during 80 d will deliver some 115 000 data points, which 
will result in a factor ~ 20 better frequency precision than for the 
Campaign 0 data. Given their scarcity and importance for stellar 
and galactic structure, we thus plan to apply for short-cadence data 
for the few available O stars in K2’s future FoVs of Campaigns 11 
and 13. 
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APPENDIX A: MASKS FOR LIGHT CURVE 
EXTRACTION 
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Figure Al. Our adopted mask (left) versus the one used by Vanderburg (right) for the light curve extraction of EPIC 202060091 as indicated in pink. The 
position of the star is indicated by the blue cross. 
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Figure A2. Same as Fig. lAll but for EPIC 202060092. 
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Figure A3. Same as Fig. lAll but for EPIC202060093. For this target, the mask on the left excluded the central saturated column to produce the photometry 
shown as black circles in the top panel of Fig.^ while this column was used to produce the Vanderburg version of the light curve shown as grey crosses in the 
top panel of Fig.^ 
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Figure A4. Same as Fig. lAll but for EPIC 202060097. 
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Figure AS. Same as Fig. lAll but for EPIC 202060098. 
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